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ABSTRACT 

Some of the most enduring questions in physics — including the quantum measurement problem and 
the quantization of gravity — involve the interaction of a quantum system with a classical environment. Two 
linearly coupled harmonic oscillators provide a simple, exactly soluble model for exploring such interaction. 
Even the ground state of a pair of identical oscillators exhibits effects on the quantum nature of one oscil- 
lator, e.g., a diminution of position uncertainty, and an increase in momentum uncertainty and uncertainty 
product, from their unperturbed values. Interaction between quantum and classical oscillators is simulated 
by constructing a quantum state with one oscillator initially in its ground state, the other in a coherent 
or Glauber state. The subsequent wave function for this state is calculated exactly, both for identical and 
distinct oscillators. The reduced probability distribution for the quantum oscillator, and its position and mo- 
mentum expectation values and uncertainties, are obtained from this wave function. The oscillator acquires 
an oscillation amplitude corresponding to a beating between the normal modes of the system; the behav- 
ior of the position and momentum uncertainties can become quite complicated. For oscillators with equal 
unperturbed frequencies, i.e., at resonance, the uncertainties exhibit a time-dependent quantum squeezing 
which can be extreme. 

PACS number(s): 03.65.-w, 04.60.-m 
I. INTRODUCTION 

Interaction of the quantum and classical worlds — i.e., between a system for which quantum effects 
predominate and a system operating in a classical limit — lies at the heart of some of the most profound and 
vexing problems in physics. The quantum measurement problem, of course, is one of the central mysteries 
of quantum physics; it has engendered concepts from Schrodinger's Cat and the many-worlds interpretation 
to consistent histories. The search for a quantum theory of gravitation also entails linking quantum matter 
sources to apparently classical spacetime geometry, as the Einstein field equations display. While now largely 
taken for granted, the notion that there must exist a quantum version of gravitation was at one time the 
subject of some debate [1], and even of experimental test [2]. One longstanding argument that gravitation 
must be a quantum phenomenon with a classical limit, and not a strictly classical field, is that an ordinary 
quantum system coupled to a truly classical field would "radiate away" its quantum nature — commutator 
values, uncertainties, etc. — leading to observable violations of quantum mechanics. A pair of simple harmonic 
oscillators with linear coupling provides an ideal test bed for exploring this aspect of quantum/classical 
interaction. 
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1 E-mail: rachmcd@mit.edu 

2 E-mail: redmouih@slu.edu 
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Two linearly coupled quantum- mechanical simple harmonic oscillators, e.g., two masses on springs, 
connected by a third spring, constitute an ideally simple model of quantum/classical coupling. First, the 
system separates into normal modes behaving as independent oscillators, so the evolution of the system from 
any initial data can be followed exactly. Second, the classical limit of a quantum oscillator is easily described 
by a coherent or Glauber state [3,4]: a Gaussian wave packet of fixed width, the centroid of which follows a 
classical trajectory. In the present work we analyze the behavior of coupled oscillators with one initially in 
its quantum ground state, the other initially in such a coherent state. The problem of a quantum oscillator 
coupled to an actual classical, i.e., c- number oscillator calls for different techniques and will be examined in 
a subsequent paper. 

The simplest version of the two-oscillator system consists of two identical oscillators, with equal 
masses, spring constants, and frequencies, plus a connecting spring with its own spring constant. Even in 
the simplest of quantum states — its ground state — this system reveals interesting effects on the oscillators' 
dynamics: The position uncertainty of an individual oscillator is reduced below its uncoupled value, while 
its momentum uncertainty is increased. (Hence, the second oscillator does not act as a "heat bath," which 
would increase both uncertainties.) The product of the uncertainties is increased; the Uncertainty Principle 
is not violated. If the system is started with one oscillator in its quantum ground state and the other in a 
quasi-classical coherent state, e.g., if the coupling between them is "turned on" at some initial time, then 
subsequently the two normal modes evolve one in a coherent state, the other in a modified state termed a 
displaced squeezed state, e.g., in quantum optics and the analysis of gravitational-wave detectors [4,5]. The 
initially quantum oscillator acquires an oscillating position expectation value — a "beat" between the normal 
modes. Its position uncertainty oscillates through values below its uncoupled, ground-state value, while its 
momentum uncertainty oscillates above its uncoupled value. The product of these uncertainties oscillates 
through values above U/2, never violating the Uncertainty Principle. 

The general two-oscillator system, with arbitrary masses, spring constants, and frequencies, can like- 
wise be analyzed exactly. We construct the wave function for a state with one oscillator initially in its 
ground state, the other in a coherent state. From this we obtain position and momentum expectation values 
and uncertainties for the first oscillator. Expectation values behave very like those in the identical-oscillator 
("symmetric") case, but the expressions we find for uncertainties in the fully general case are rather opaque. 
However, one simple case may be of particular physical importance: that of oscillators with different masses, 
but equal (uncoupled) frequencies, i.e., that of quantum and classical oscillators interacting "at resonance." 
In that case quantum uncertainties clearly exhibit the behaviors found in the symmetric case: Position 
uncertainty is reduced from its uncoupled quantum value-toward the classical zero value-while momentum 
uncertainty is increased. The product of these remains always above the Uncertainty Principle bound; the 
system of course never actually violates quantum mechanics. 

The simplest case of two identical oscillators is detailed in Sec. II below. The fully general case is shown 
in Sec. Ill; the special but physically important equal-frequency case is treated in Sec. IV. We summarize 
our results and conclusions in Sec. V. For ease of comparison between quantum and classical regimes, we 
retain the constant h explicitly throughout. 

II. SYMMETRIC COUPLED OSCILLATORS 

The simplest system of two coupled harmonic oscillators consists of two identical, one-dimensional 
oscillators with linear coupling. For example, two identical masses m, each on a linear spring with spring 
constant k [so each oscillator has uncoupled angular frequency u> = (fc/m) 1 / 2 ], connected by a third spring 
with spring constant k, constitute such a system. 
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A. Hamiltonian and normal modes. 



The Hamiltonian for this system can be written 



H=^^ + \kx\ + \kxl + \K{x 2 -x 1 f (2.1a) 



= P+ ^~ + h kx + + h( k + 2k ) x - > ( 2 - lb ) 

in terms of individual oscillator displacements (from equilibrium) x\ and x 2 and momenta p\ and p 2 , or 
normal- mode coordinates x± and momenta p± . The normal-mode amplitudes are related to the individual 
coordinates via 

x + = ^ (2.2a) 

- ^ , (Mb) 

where the scale factor is chosen to make this a unitary transformation. The normal-mode momenta p±, of 
course, are conjugate to these. The separated form (2.1b) of the Hamiltonian allows the system to be treated 
as two independent harmonic oscillators. 

B. Ground state. 

The quantum ground-state wave function of the system can be written as the product of the normal- 
mode ground states, viz., 

. . /m^fA 1 ^ 4 / mui 9 \ / mfl , \ r ,• . ^. n . . 

*o(x±,t) = y-^Tj exp(-— 4J cxp^-— z 2 _J cxp[-f(w + f!)i] . (2.3) 

The normal- mode angular frequencies follow from Hamiltonian (2.1b): 

/ k \ 1/2 

w = - (2.4a) 



for the + mode, and 

n/t+M'".-, (2 . 4b) 

V m J 7 
for the — mode, where the last expression defines the parameter 7. 

The probability distribution for the position of one oscillator, with the system in this state, is obtained 
by integrating the probability density over the other oscillator coordinate. The result is 



/oo 
%{x ± ,t) V (x±,t) dx 2 
-CO 

fm 2 Lon\ 1/2 r ( 

= U^J L cxp { 



-^-(zi+zs) -— (z 2 -a;i) ) cfe 2 



2mw \ / 2muj 



nh(l+^)J Ml +7) 
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e *P I s *i) ■ ( 2 - 5 ) 



This is a Gaussian probability distribution with standard deviation 



<7« 



2mu) 



^ • (2.6) 



The unperturbed vacuum-state probability distribution for the position a; of a single oscillator is 

^t) = {—) ex P (-— (2.7a) 

with standard deviation 



With ft > uj, i.e., 7 < 1, the effect of the coupling is to reduce the position uncertainty of the oscillator. 

The probability distribution for the momentum of one oscillator can be calculated similarly or by 
utilizing the symmetry of a harmonic oscillator under the interchange y/rruv x <-» pj y/muj. The result is 

1 /2 

Pi(Pi.*)= ( ^J 1 ! , expf- ^— ^ ; - p{] , (2.8a) 

with standard deviation 



\irhmuj(l + 7) / \ frmu(l + j) 



These can be compared with the unperturbed results 

1 \ V2 



<t<°> 



hmuj 



(2.9b) 
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The coupling increases the momentum uncertainty of the individual oscillator. 

Thus the net effects of the coupling, in the ground state of the coupled system, are: 

• To decrease the position uncertainty of the individual oscillator; 

• To increase the momentum uncertainty of the oscillator; 

• To increase its uncertainty product a x a p above the minimum value h/2; 

• To break the symmetry between y/mujx and p/yjmu for the individual oscillator. 

The last is not so unexpected, since the coupling perturbation is in position and not symmetrically in 
momentum. 

The effects of the coupling cannot be characterized by the introduction of a finite temperature for the 
individual oscillator. That is, the second oscillator does not act as a "heat bath," even though coupling to 
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a large number of oscillators should do so. The probability distributions and uncertainties for an oscillator 
in thermal equilibrium at temperature T are [6] : 



V t(x) = tanh(/3/2)) 1/2 cxp (-^ tanh(/3/2) x 2 ) 



— coth (/3 /2) 



^ tanh(/3/2) V /2 / tanh(/3/2) 2 N 

7>t(p = r exp p z 

\ irnmuj J \ nmui j 



a, 



(T) _ 



fimuo 



coth(/3/2) , 



(2.10a) 
(2.10b) 
(2.10c) 

(2.10d) 
as well as the 



with j3 = hio/(kBT). The net effects of a finite temperature are to increase both a x and a 
uncertainty product <j x <j Pi while preserving the symmetry between y'wx and p/y/mu). 

C. Coupled ground/coherent state. 

With the system in its ground state, both oscillators display fully quantum dynamics. We model a 
quantum oscillator coupled to a classical oscillator by a different choice of state: We envision oscillator #1 in 
its unperturbed ground state, and oscillator #2 in a coherent or Glauber state with classical amplitude X , 
at time t = 0. We follow the subsequent evolution of this state, and determine the probability distributions, 
expectation values, and uncertainties for oscillator #1. In effect, oscillator #1 is in its ground state, oscil- 
lator #2 is behaving classically with oscillation amplitude X n , and the coupling "turns on" at time t = 0. 
For simplicity, we assume the classical momentum of the second oscillator is zero at t = 0, i.e., it is at a 
maximum of its classical trajectory at that instant. 

1. Initial wave function. 

The initial wave function for the system, i.e., the wave function at time t = 0, for this coupled 
ground/coherent state is [3,4]: 



, /muj\ L ' z j muj 2 \ / rruv . 

nxuX2 , 0) = j exp (-— Xl ) exp (- — (x 2 X y 

In terms of the normal-mode amplitudes, this is simply 



(2.11a) 



/"muj\ 1 / 2 

* (x± ' 0) = br) exp 



7TILO I Xq 



2h 



cxp 



2h 



V2 



(2.11b) 



This corresponds to a coherent state of the + mode, with initial amplitude Xq/^/2 and momentum zero. The 
initial state of the — mode, however, is a displaced squeezed state [4,5]: The initial displacement is Xo/V2, 
the initial momentum is zero, but the initial width of the Gaussian wave packet is proportional to l/\/u> 
rather than the 1/VU required for a coherent state of this mode. (In fact this is a displaced "antisqueezed" 
state, with l/y/tu > 1/VU.) 
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2. Propagation of the wave function 



The wave function for all times t > is obtained by applying harmonic-oscillator propagators [7] for 
the normal modes to the initial wave function ^(x± } 0). The result is the double integral 



/OO pOC 
J 
-co J — 



sin(wt) 



1/2 



-co J —CO 



x exp 



x exp 



mio 
IE 

imuj 



27n?isin(£lt) 

2 



X 



1/2 



exp 



Xo 
V2 



2hsm(tut) 
2hsiri(nt) 



[{y 2 + + x 2 + ) cos(wt) - 2x +y+ ] J 
[(y 2 . + x 2 _) cos(fii) — 2x-yJ\ I 



(2.12) 

where coordinates y, at arbitrary time t, are distinguished from coordinates x at t = 0. The integral over x + 
gives a coherent state in y + , with classical amplitude (Xo/y/2) cos(wi) and momentum -mw(Xo/\/2) sin(wi). 
The integral over x_ gives a displaced squeezed state in y_. Ultimately, the wave function takes the form 



*(y±,i) = ipc(y+,t)ipds(y-t) , 



(2.13a) 



with [3,4] 



/mu)\ 



1/4 



exp 



exp 



47i 



■Xq sin(wt) cos(u>t) 



x exp 



imw Xq 



sin(u;i) y + I exp 



-2H [ y+ - V^ C ° S{UJ ) 



(2.13b) 



and [4,5] 



tpds(y- 



/mwx 1 / 4 / cos(Qt) — ijsiriiilt) \ ^ 2 ( imujXQ cos(ftt) sin(ftt) \ 
. f ) = J I —777^ , — 27777; I ex P I 777— ?77^ , ~ 27777^ I 



x exp 



x exp 



cos 2 (m) + 7 2 sin 2 (m) 

-imw7(X /v / 2) sin(fit) y. 
h[cos 2 (nt) + 7 2 sin 2 (m)] 



exp 



y- 



^0 

V2 



4^[cos 2 (m) +7 2 sin 2 (m)]y 

«mO(l - 7 2 ) cos(fii) sin(Qi) y 2 
2h[cos 2 (Qt) + 7 2 sin 2 (fit)] 

21 



cos(fii) 



(2.13c) 



2^[cos 2 (m) +7 2 sin 2 (m)] 

in terms of normal-mode coordinates y± and time t > 0. 

5. Reduced probability distribution for y±. 

As before, the probability distribution for the single oscillator coordinate y\ is obtained by integrating 
the full probability density over y 2 - We obtain 



/oo 
\ip c (y+,t) tp ds (y-,t)\ 2 dy 2 
-OO 



6 



nh V cos 2 (fit) + 7 2 sin 2 (m) 



1/2 



/ 



x / exp 



muj { 



exp 



( ( x ° ro+y\ 



h cos 2 (fit) + j 2 sin 2 (fit) 



V 



J 



2 ma 



7r7i[l + cos 2 (m) + 7 2 sin 2 (m)] 



1/2 



exp - 



2ma; {yi - ^X a [cos(Lut) - cos(fit)]} 2 
H 1 + cos 2 (fit) +7 2 sin 2 (fit) 



2 ma; 



7rfi[l + cos 2 (fit) + 7 2 sin 2 (fit)] 



1/2 



exp 



2 ma; 



Sl-u \ / fi + a; 
2/i — Xq sin I — - — t ] sin | — - — t 



1 + cos 2 (fit) +7 2 sin 2 (fit) 



(2.14) 

treating y + and y_ as functions of forms (2.2a) and (2.2b) of coordinates y\ and y 2 - 

4- Quantum expectation values and uncertainties for oscillator #1. 

The expectation value and uncertainty for the position of oscillator #1 can be read directly from the 
Gaussian distribution V(yi,t). We find 



. . X , , . . / fi — u) \ . ( fi + u 

\2/i) = — [cos(u;i) - cos(fii)J = X sin I — - — t I sin I — - — t 



and 



r(!) 



4mw 



[l + cos 2 (fit)+ 7 2 sin 2 (fit)] . 



(2.15a) 



(2.15b) 



The behavior of (yi) can be quite complicated; it can be as large as Xq or as negative as — Xq. In the case 
of weak coupling, i.e., k <C k or fi — u <C to, it oscillates with a slowly oscillating amplitude, exhibiting 



"beats" between the normal modes. The position uncertainty off oscillates between the unperturbed value 
y/h/ (2muo) and the smaller value \/h(l + 7 2 )/(4ma;). 

Momentum expectation values and uncertainties cannot be determined from the reduced probability 
distribution V(yi,t); the full wave function ^(y±,t) is needed. The normal-mode relations (2.2a) and (2.2b), 
plus the chain rule, imply 



Pi = 



P+ -P- 



whence follow 



and 



(Pi) 



V2 ' 
~ (P-) 

L(+) 2 , r-)2 



(2.16a) 



(2.16b) 



(2.16c) 
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the last since the normal modes are uncorrelated in this state. The necessary expectation values and 
uncertainties can be read from the momentum-space wave functions for the normal modes, viz., 



■>PdP+,t) 



oo e -i p+y+ /h 



i 



mivirh 



1/4 



exp 



exp 



imui 



Xfi cos(wt) sin(wt) 



i -X.Q . 
x exp [ — J^~^ cos(u)t)p+ J exp 



1 



2mujh 



p+ + mu—^L sin(o>t) 
V2 



and 



-00 



,(?-,*) = J 



00 e -ip-y-/K 

(2^)V2 



ipds{y-,t)dy- 



_/ 7 y/4 / 7 cos(fit) - isin(fit) \ 1/2 / -imtl^Xl cos(fit) sin(fif) 
" ImfiTrftJ V7 2 cos 2 (fit) + sin 2 (fit) J V 47i[ 7 2 cos 2 (fit) + sin 2 (fit)] 



x exp 



-17 2 (X /V2) cos(fit)p_ 
ft[ 7 2 cos 2 (fit) + sin 2 (fit)] 



exp 



/ i(l- 7 2 )cos(fit) &m{Q,t)p 2 _ \ 
\2mnh\~f 2 cos 2 (fit) + sin 2 (fit)] ) 



(2.17a) 



These imply 



x exp 



2mfi^[ 7 2 cos 2 (fit) + sin^(fit)] 



p_ + mfi— p sin(fit) 
V 2 



(2.17b) 



and er^ ' 



(p+) = 


-mw-p sin(wt) , 
v2 


(2.18a) 


„(+) = 

p 


/ mwli 
V 2 ' 


(2.18b) 


(p-> = 


— mfi^^ sin(fit) , 
v2 


(2.18c) 



[ 7 2 cos 2 (fit) + sin 2 (fit)] 



(2.18d) 



Hence, relations (2.16b) and (2.16c) yield the results 

ttlXq 



and 



<Pi> 



<7« = 



cj sin(wt) — fi sin(fit)] 



' mu>% 



1 + cos 2 (fit) + — sin 2 (fit) 



(2.19a) 
(2.19b) 



The momentum and position expectation values are related "classically," via (pi) — m d(y\) jdt. The momen- 
tum uncertainty oscillates between the unperturbed value yJmojfi/2 and the larger value ^muh(l + 7 _2 )/4. 



Thus the behavior of oscillator #1 is perhaps reminiscent of the suggestion that a quantum oscillator 
coupled to a classical one would lose its quantum nature. In particular, its position uncertainty is reduced 
below its unperturbed ground-state value. It cannot, however, be reduced arbitrarily close to the classical 
zero value. For any value of the coupling, i.e., of 7 — and independent of the classical amplitude Xq — the 
uncertainty is never less than l/\/2 of the unperturbed quantum value. The effect of the coupling could 
more aptly be described as a quantum "squeezing" than a loss of quantum nature. (Of course, since this 
model remains a quantum-mechanical system no matter what its state, it could never exhibit an actual 
violation of quantum mechanics.) The uncertainty product for oscillator #1 is 



v p 




cos 2 (fit) 



■ 7 2 sin 2 



(fit)] \l 



1 



cos z (fit) + — sin^(fit) 



U + (I_ 7 )1^!M 



m7 +7 ) ~G~ 7 ) cos4(!!<) - 



This oscillates between the values 



!<4'H"<!(i +7 



(2.20) 
(2.21) 



always in accord with the Uncertainty Principle. 

III. GENERAL COUPLED OSCILLATORS 

The preceding case of two identical oscillators is rather special. One might expect coupled quantum and 
classical oscillators to have very different masses and spring constants — including, for example, a quantum 
oscillator coupled to a mode of a radiation field. The above analysis is readily generalized to the case of 
oscillators with arbitrary masses mi and m 2 and individual spring constants k x and k 2 , hence unperturbed 
angular frequencies u\ = \Jk\jm\ and ui 2 = \fk 2 /rn 2 . 

A. Normal modes. 



The Hamiltonian for this two-oscillator system, coupled as before by a spring with spring constant k, 



is 



H 



1 + -^ 2 — + hkixj + \k 2 x\ + \k{x 2 - xi) 2 . 



2mi ' 2m 2 ' 2 1 1 2 ' 2 2 
To separate this into normal modes, we first rescale the coordinates and momenta thus: 



The Hamiltonian becomes 



X x = 
Pi = 

H = 



mi 
m 2 

m 2 
mi 



1/4 



Xi 



1/4 



Pi 



x 2 =[^ 

mi 



m 2 



1/4 



x 2 



1/4 



P2 



(3.1) 

(3.2a) 
(3.2b) 



+ 3« 



\ 1/4 / \ 1/4 

mi \ ^ I m 2 N 
m 2 J 2 \mi 



X 1 



(3.3a) 



with geometric-mean mass 



fi = (mim 2 ) 1 / 2 

Now the rotation of coordinates and momenta 



hence, 

transforms the Hamiltonian into 



x+ 

X- 

P+ 
P- 



cos a sin a 
— sin a cos a 



cos a sin a 
- sin a cos a 



X 1 

x 2 

Pi 
p 2 



H 



given the conditions 



a = 2 arctan 



2k/ h 



ui-uf 1 - 



(3.3b) 

(3.4a) 
(3.4b) 

(3.5) 
(3.6a) 



UJ+ 



2 2 2*2 

w, cos a + Wo sin cH — 



mi 
m 2 



1/4 



sin a 



m 2 
mi 



1/4 



-i 2 



cos a 



1/2 



(3.6b) 



and 



sin 2 a + mi, cos 2 a 



1/4 



/ mi \ ' / m 2 

cos a + — | sin a 

\m 2 J \nii 



1/4 



1/2 



(3.6c) 

fixing the rotation angle a and the normal- mode frequencies u±. With this Hamiltonian, the normal 
modes x± evolve as independent harmonic oscillators. The combined transformations (3.2a) and (3.4a), 
between (x\ 1 x 2 ) and (x + ,X-), remain a unitary transformation in the quantum- mechanical sense, though 
not as a matrix for mi ^ m 2 . 

B. Coupled ground/coherent state. 

As above, we model a quantum/classical coupling via a state in which oscillator #1 is in its ground 
state, and oscillator #2 is in a coherent state, with classical amplitude x and momentum zero, at initial 
time t = 0. The initial wave function is 

1/4 



/mim 2 wiw 2 \ / m 1 uj 1 xl+m 2 uj 2 (x 2 - x ) 2 \ 

*&>**>°) = { -itf ) ex H m J ' 



(3.7a) 



in terms of the original oscillator coordinates, or 



/ 2 \ V 4 / 

-r / ^\ / M ^1^2 \ / M / 9 • 9 \ 9 / • 9 

w(x±, 0) = 2— I exp — — (ivi cos a + W2 sin a)x^_ + (wi sin a + < 

V 7r 2 7i / \ 2n L 



UJ 2 cos 2 OL)X 2 _ 

+ (uj 2 — uji) sin(2a)x + a;_ — 2w 2 A (sin a x + + cosaa;_) + uj 2 Xq 



(3.7b) 
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in terms of normal-mode coordinates, with X = (rr^/mi) 1 / 4 ^ a rescaled classical amplitude. 



The time-dependent wave function for this state is again obtained by applying separate harmonic- 
oscillator propagators G± for the two normal modes. The integrand of the propagation double integral takes 
the form 



G+G_*(x ± ,0) 



2 \ V 4 



-k 2 H 2 ) \47r 2 i 2 ?i 2 sin(w + i) sin(o)_t) 



1/2 



exp 



2h 



x exp 



(ifiw + cot(u + t) 2 ^ ( iflUJ- cot(w_t) 

2h~ 



V+ exp 



2h 



x exp | loi cos 2 a + lu 2 sin 2 a — iuj + cot(u> + t) x 2 + 

+ loi sin 2 a + w 2 cos 2 a — iuj- cot(u_t) x 2 _ + (lu 2 — loi) sin(2a) x + x_ 

— 2 LU2X0 sin a — iu + csc(ui + t) y + x + — 2 o>2^o cos a — csc(a>_t) y_J X- 

(3.8) 

Hence, the propagation integral over x± is a two-dimensional Gaussian integral of the form 

j cxp (_|_(A T <M - D T A - A T D)) d 2 x = J exp (-^(B T SB - D T S- 1 D)) d 2 ^ 



4n 2 h 2 N 1/2 



det S 



with 



W2^o sin a — iuj + csc(ui+t) y+ 



D = 



and 



U2X0 cos a — iuo- csc(w_i) y_ 



' lui cos 2 a + w 2 sin 2 a — iuj + cot(uj + t) 
^(lo 2 - wi) sin(2a) 



(3.9a) 

(3.9b) 
(3.9c) 

(3.9d) 



\{ui2 - lui) sin(2a) 
u>i sin 2 a + LU 2 cos 2 a — iuj- cot(w_t) 



(3.9e) 
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The wave function is then 
*(»±,t) = 

with 
and 



2 \ V4 

(1 LO\UJl \ 



Tr 2 h 2 J \i 2 sm(u)+t) sin(w_t) det S 



1/2 



exp 



M^2^Q 

27i 



x exp + . y r 02 



JJ+ 

y- 



( LJ-L. cot(w+i) 





o>_ cot(w_i) 



for all i > 0. Rewriting column matrix (3.9d) as 

D = -iSl^Y + iZ) 

with 



o — ( UJ + csc(ai + £) 
" 1= ^ 





w_ csc(w_t) 



and 



Z = 



puts the wave function into the form 



1/4 



^2 

— X sina sin(cj + i) 
— Xq cos a sin(u;_i) 



i 2 sin(w+t) sin(w_t) detS 



1/2 



exp 



M w 2^0 



x exp (-^ \(Y T + iZT^S^VL^Y + iZ) - iY T fl 2 Y] ) 



in terms of the matrices defined above. 



(3.10a) 



(3.10b) 



(3.10c) 



(3.11a) 



(3.11b) 



(3.11c) 



(3.12) 



To determine probability distributions, expectation values, and uncertainties, it is useful to separate 
the argument of the exponential in $ into real and imaginary terms. The first matrix of coefficients can be 
written 

= U + iV , (3.13) 
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where U and V are real, symmetric 2x2 matrices, with U nonsingular and positive definite. The wave 
function is then 



i 2 sm(uj + t) sin(w_i) det S 



1/2 



exp 



2h 



XeXP "2l 



Y T UY - Z T VY - Y T VZ - Z T UZ 



+ iY T (V - tt 2 )Y + i(Z T UY + Y T UZ) - iZ T VZ 



2 \ V4 



z 2 sin(w + i) sin(w_i) det S 



1/2 



exp 



2h 



X8Xp "2* 



(Y T - Z T VU- 1 )U(Y - U^VZ) - Z T VU- 1 VZ - Z T UZ 



+ iY T (V - Q 2 )Y + i(Z T UY + Y T UZ) - iZ T VZ 



(3.14) 

This can be simplified: Since the matrix U contains no dependence on the classical amplitude Xq, normal- 
ization of the wave function implies 



oj 2 X 2 = Z T (VU- 1 V + U)Z , 



and also 



det 17 = 



U\U 2 io\u) 2 _ 



sm 2 (LU+t) sm 2 (uj-t)\ detS\ 2 



(3.15a) 



(3.15b) 



These can also be verified using the explicit expressions shown in Sec. D below. The first result reduces the 
wave function to the form 



*(y,t) 



2 \ V 4 



1/2 



n 2 h 2 J \i 2 sm(uj + t) sin(w_i) det S 

x exp [Z T VZ - (Z T UY + Y T UZ) - Y T (V - Ct 2 )Y 

xexp(-^(r T - Z T VU- 1 )U{Y -U^VZ)) . 



(3.16) 



This wave function is finite or regular for all y± and t > 0; the components of U, U^ 1 , V — Cl 2 , and VZ 
contain no vanishing denominators or divergent circular functions. This too can be verified explicitly using 
the expressions in Sec. D below. 

The probability distribution for the normal-mode coordinates is 
P(Y,t) = \*(Y,t)\ 2 



V 2 det U 
n 2 h 2 



1/2 



exp (~j-(Y - U-^VZfUfy - U^Vzfj 



(3.17) 
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a Gaussian distribution in two dimensions. To obtain the reduced probability distribution for the single 
oscillator coordinate y±, the inverse of the transformation to normal- mode coordinates, 



Y 



cosa sina^ f (toi/to 2 ) 1 ^ 4 




sin a cos a 







(to 2 /toi) 1 ^ 4 



Ny 



(3.18) 



is needed. Here y denotes the column containing the oscillator coordinates, and N the product of the two 
square matrices. The matrix N has unit determinant, but is not orthogonal if masses m\ and m 2 are unequal. 
Distribution (3.17) can be expressed thus: 



p(y,t) = 



fi 2 det U 



1/2 



cxp (-|CV - N^U- 1 VZ) T N T UN{y - N^U^VZ^j 



(3.19) 



as a distribution for y\ and y 2 . Reduction to V(yi,t) requires an integral of the general form 
J cxp — (xi x 2 ) ^ (jj, 1 ^ d x i= J cxp(— ax\ — 2cx\Xi — bx 2 ) dx 2 



(?) 



1/2 



exp 



ab - c 2 



(3.20) 



Hence, integrating V(y,i) over y 2 yields 

/zdetZJ ^ 1/2 



irh(N T UN)r 



exp 



/i det U 



h(N T UN) 2 



yi - (JV- 1 ^- 1 ^)! 



(3.21) 



a normalized, Gaussian distribution for y\. 

C. Expectation values and uncertainties. 

The expectation value and uncertainty for the position of oscillator #1 in this state can be read off 
of V(yi,t). They are 



( yi ) = (N-'U-'VZ), 



and 



in terms of the above-defined matrices. 



h (N T UN) 22 



2fi det U 



(3.22a) 
(3.22b) 



To calculate the momentum expectation value and uncertainty, the full wave function (3.16) must be 
used. First, we obtain 



(Pi) = / ^*-^-^d yi dy 2 
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T 



J |( 1 )N T UN(y - N^U^VZ) + {y- N^U^VZf^UN (^j 



+ 1 



Z T UN ( g ) + (1 0)N T UZ 



(i o)Ar J '(y-fi 2 )ivy + y T iv T (y-o 2 )Af 



P(y,t)d yi dy 2 



= -^{n t [u + (v-n 2 )u- 1 v]z} 1 , 

this last following from the matrix expectation value (y) = N~ 1 U~ 1 VZ . Second, we find 



(3.23) 



0)7V T lw(M+z(l 0)iV T (T/-f7 2 )ivP 



+ ^2 ({ [^ T ^(^ - A^l^VZ)] x + i [N T {V - fi 2 W] 1 + li^UZ^}' 

= ^h{N T [u + i(v-n 2 )}N} 11 
-n 2 (^(y- n- x u- 1 vz) t n t [u + %{v - n 2 )]N ° ^n t [u + i(v- n 2 )]N(y - n^u^vz) 

+ ,i 2 ({N T [U+(V- Q^U-'V] Z}^ 2 , (3.24) 
using the facts the square of the first component of any column matrix C can always be written 



c 



(3.25) 



and that terms linear in y — N 1 U X VZ appearing in the original squared expression have zero expectation 
value. The last term in Eq. (3.24) is precisely (pi) 2 — without specifying the matrices — leaving 

4 1)2 - (pi) (pi) 2 

= fih{N T [u + i{v -n 2 )]N} n 

-v 2 ((y-(y)) T N T {u + l (v-n 2 )}N( 1 Q °) n t [u + t (v -n 2 )]N(y - (y))^ . 

(3.26) 

The last expectation value is the integral of a quadratic form times the normalized Gaussian V(y,t). The 
general form of such an integral is 



^^y /2 | X T AXe- xTBX d n x= ^{AB- 1 ) , 



(3.27) 



15 



for positive-definite B, as can be proved by diagonalizing B or by explicit integration of component expres- 
sions. This implies 

= fih{N T [u + i(v - n 2 )}N} n 

N T [U + i(V-fl 2 )}N ° ^jN T [U + i(V-n 2 )}N^N T UN^ 



M 2 T 
~2 



= y i NT l U +{ y- ^) U ~\ V - fi 2p} n . (3.28) 

As expected, this is real and of order Ti, even without using the specific forms of U, N, V, or £1 2 - It implies 
the momentum uncertainty 



(1) /M 



{N T [U +(V- n 2 )U-HV - fl 2 )]N} u . (3.29) 

The probability distribution for momentum pi is a Gaussian, with expectation value (pi) and standard 
deviation . 

D. Explicit expressions. 

Explicit expressions for the wave function, probability distributions, expectation values, and uncer- 
tainties in the coupled ground/coherent state — in terms of the individual oscillator masses, frequencies, and 
coupling — involve an unwieldy number of terms, even with the aid of a computer. It facilitates matters to 
define the following real and imaginary parts, in terms of the matrices in Eqs. (3.9e), (3.10c), and (3.11b): 



with real matrices 



and 



and frequency combinations 



m = n^sn^ 1 = n- a , (3.30a) 



uo c p 2 8 sin a cos a prj N 

Tl=\ | (3.30b) 

5 sin a cos a prj u> s r] 2 



' pcos(uj + t) 

I=Q^ 1 n 2 Q^ 1 = ( ) , (3.30c) 

rjcos(uj-t) 



uj c = lji cos 2 a + uj 2 sin 2 a , (3.30d) 

uu s = u>i sin 2 a + uu 2 cos 2 a , (3.30e) 

S = uj 2 — coi , (3.30f) 

p = , (3.30g) 

sin(aj_f) 

and r] = — . (3.30h) 



The inverse matrix M 1 is the matrix in Eq. (3.13), implying 

1 

M 



U + iV = -J- (K - il) , (3.31) 



1G 



with M. = det M and tildes denoting adjoint matrices (transposed matrices of cofactors, not Hcrmitian 
conjugates). The matrices appearing in the wave function, et cetera, can then be written: 

U = \J^ mM)7Z { ^ M)i] ' (3 - 32a) 

U ~' = \MftetU mM)n { ^ M)I] ' (3 - 32b) 
V = -^[(3A^)K + • (3.32c) 

The real and imaginary parts of M. are, explicitly, 

= prj[uiuj2prj — cos(ui + t) cos(w_t)] , (3.33a) 
and — — pt][u) s r] cos(cu + t) + uj c pcos(ui^t)] . (3.33b) 

Identity (3.15b) is equivalent to 

\M\ 2 dct U = u} 1 Lu 2 p 2 r l 2 =dctK , (3.34) 
which can be confirmed by direct evaluation. Finally, 

(n-iZ)(n-ii) = Ml , (3.35a) 

where 1 denotes the 2x2 identity matrix, implies the identities 

$M = det "K- det I , (3.35b) 

as may also be obtained directly, and 

1K + K1= -SM 1 . (3.35c) 
These simplify the calculations considerably. 

Evaluation of the position expectation value (yi) is straightforward. Results (3.32b) and (3.32c) and 
identities (3.35b) and (3.35c) imply 

U ~ 1V = jMjdrtU 171 ■ (3 ' 36) 

Explicit evaluation yields 

' r\ sin a 

KZ = uj^Xopr] I I . (3.37) 



pcosa 



With identity (3.34) and form (3.30c), this gives 



Xo sin a cos(uj + t) 

(Y) = XJ- X VZ = ( ) , (3.38a) 

Xocosa cos(w_t) 



whence follows 



(y) = n- x u- x vz 

(rri2/mi) 1 / 4 X sin a cos a [cos(uj + t) — cos(w_i)] 
(mi/m 2 ) 1 / 4 X [sin 2 a cos(w+i) + cos 2 a cos(w_t)] 



(3.38b) 
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This gives the desired result: 

(yi) — { m 2/ m i) 1 ^ 2x Q sin a cosa [cos(w + t) — cos(o>_i)] , 



(3.39) 



in terms of the original classical amplitude x$ and angle a and normal-mode frequencies oj± from Eqs. (3.6a- 
c). 



Evaluation of the position uncertainty requires the matrix 

N T UN = ——-[(31M)N T 1ZN— (^sM)N T TN] . 



\M\ 



Direct evaluation of the matrix elements 



1/2 



(N J RN) 22 = ( — - ) (uJsV 2 sm2 a ~ ^Sprj sin 2 a cos 2 a + uj c p 2 cos 2 a) 
\rnxj 



and 



rn 2 



1/2 



(N TN)22 = [77 sin a cos(o>_i) + pcos a cos(u; + i)] 

\m 1 > 



enables us to obtain the final result 



cos a — cos {ui + tj H 5— sin (w+i) 



U! 2 



1 - ■ 2, mi w l u s .2/ 

+ sin a — cos (w_£) H *— sin (w_i) 



^2 



On 9 / ^1^2 

+ 2 — cos a sin a I cos(u;+f) cos(w_t) sin(a; + i) sin(w_i) 

W2 \ LO+U)- 



1/2 



This rather complicated result reduces to form (2.15b) in the equal-mass, equal- frequency case. 



The momentum expectation value (pi) is obtained from the matrix 

1 

\M\> 



U + (V- n 2 )U- l V = [^M)TZ (ZM)i\ 



1 



\M\ 



:[(9fM)^+(»M)2]+fi: 



\M\ 2 &etU 



in 



\M\ 2 AetU 



WlLU2p 2 r] 2 



(1 - n 2 i)n 

sin 2 (w + t) 



sin 2 (w_i) 



Result (3.37) then yields the column matrix 

[U + (V - o 2 )*7~V]z = X Q 



uj + sin(w + t) sin a 
w_ sin(w_i) cosa 



(3.40) 



(3.41a) 



(3.41b) 



(3.42) 



(3.43) 



(3.44) 
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hence. 



N T [U +(V - U2)U-*V\Z = X 



(mi/TO 2 ) 1 / 4 sina cosa[uj + sin(uj + t) — w_ sm(w_t)] 
_ (m 2 /'™i) 1 ^ 4 [sin 2 acj + sin(oj + i) + cos 2 sin(w_i)] 



(3.45) 



(3.46) 



Result (3.23) takes the form 

/ \ 1/2 

(pi) = —mix ( — ) sin a cos a [ui + sin(a; + i) — w_ sin(u;_i)] . 

As in the symmetric case, the expectation values (t/i) and ($>\) are related by (pi) = mi d(yi)/dt. 

The momentum uncertainty Cp 1 ^ is given by Eq. (3.29). After some manipulation, the requisite matrix 
takes the form 



u + (v - n 2 )u- 1 (v - n 2 ) = -r-^P - - rczo 2 

+ (3fcM)fi 2 ftfi2 - (3?M)fi 2 X0 2 ] . 
Hence, the uncertainty is obtained from the n component of 

N T [U +(v- si 2 )u-\v - fl 2 )]N = 



(3.47) 



dot ft 



N J 



T f sin 2 (cj+t) 

sin 2 (uj-t) 



TIN 



- N T K ( cos2 ( w +*) 



cos 2 (w_t) 



N 



+ {^M)N T n 2 KQ, 2 N - (^M)N T n 2 in 2 N 



(3.48) 



The components of the individual matrix products on the right-hand side can be obtained by explicit calcu- 
lation. The results are quite long and unwieldy, even after some trigonometric simplifications. When they 
are combined, some cancellations and further trigonometric combinations leave the final result 



r(l) 



cos a [ — cos (uj+t) H * — sin (ui+t) 

1 LOi W{W 2 



+ sin a I — cos (ui-t) H * — sin (uj-t) 



— 2 — cos asm a cos(u + i) cos(u_t) sm{u) + t) sin(w_i) 



1/2 



(3.49) 



This rather complex result agrees with Eq. (2.19b) in the symmetric case. 



Similar matrix manipulations can be used to confirm identity (3.15a) explicitly. They can also be used 
to prove the finiteness of the matrices U, VZ, V — fl 2 , and U~ x appearing in the wave function 'J', thus: 
The combinations UZ + iVZ and U + i(V — Q 2 ) can be written as finite complex matrices divided by the 
complex number E(t) = M(t)/[p(t)r)(t)]. Hence, finiteness of all components of UZ, VZ, U, and V — £1 2 
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is assured provided S neither vanishes nor approaches arbitrarily close to zero at any time. This number is 
obtained from Eqs. (3.33a) and (3.33b): 

$tE(t) = — cos(uj + t) cos(o)_i) + — - sin(u) + t) sin(w_i) (3.50a) 

UJ+UJ— 

3E(i) = — — cos(u> + t) sin(w-i) — — cos(o>_t) sin(u + i) . (3.50b) 
Consequently, S satisfies the identity 



2 2 

/(x*, *2 //^i \2 2/ \ 2/ \ Ci^^(^3^i^2 -2/ \ • 2/ \ 

w c Ws(5fi^) +witJ2(^^) =w c WsCOS (w+i) cos (cu-t) H ^ — 2 — sm sm 



CJlCJ2tJ o/ \ • 2/ \ W1W2W 2/ \ • 2/ \ 

H j — cos sm ( w -*) H 2 — cos sm • 

0J_ UJi 



(3.51) 



With 



£ = max{w c w s ,wi[ii2} (3.52a) 
and A e mm w c w s , ^ — ^ — , ^ — , ^ — f > (3.52b) 



identity (3.51) implies £ |S| 2 > A, i.e., 

|S| 2 > ^ > . (3.53) 

As required, S is bounded away from zero by a time-independent bound. Likewise, matrix U' 1 can be 
written in terms of finite matrices times S. But with 



C' = min{w c o; s , wi^} = ^1^2 (3.54a) 

, w / w c w s w 2 w| W^a^ 2 WiW 2 W 2 l 
and A =max^w c w s , 2 — ^ — j 2 — ' 2 — r ' (3.54b) 



identity (3.51) also implies (' |^| 2 < A', i.e., 

A' 



S| 2 < ^ < 00 . (3.55) 



With 5 bounded away from infinity, the finiteness of U 1 is also assured. 



The expectation values (yi) and (pi) for quantum oscillator #1 show the same beating between the 
normal modes in the general case as in the symmetric case. But the uncertainties <JjP and (Tp~\ which 
incorporate the effects of the coupling on the quantum character of the oscillator, are too complicated in the 
general case to describe easily. It is more illuminating to focus on another special case — one of particular 
physical significance. The effects of one oscillator on another will certainly be most pronounced when the 
two are at or near resonance. It will therefore be most useful to focus on the case of two distinct oscillators, 
with unequal masses mi 7^ m 2 and correspondingly different spring constants, but with equal (unperturbed) 
frequencies u)i — UJ2- 
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IV. UNEQUAL-MASS OSCILLATORS AT RESONANCE. 



The behavior of coupled oscillators of distinct masses but equal unperturbed frequencies is easily 
obtained from the general results of Sec. Ill above. In this case the oscillator frequencies and frequency 
combinations reduce to u>± = w 2 = uj c = lo s = uj and 6 = 0. 



A. Normal modes. 

The normal-mode rotation angle a for this case follows from Eq. (3.6a), i.e., 



tan(2a) = — — , (4.1) 



which implies 

cos(2a) 
sin(2o:) 

hence 

cos a = 
sin a = 



mi — m 2 



mi — m 2 
mi + m 2 

2/i 
mi + m 2 



mi 
mi + m 2 



m 2 
mi + m 2 



(4.2) 



(4.3) 



These reduce to a = tt/4 in the symmetric case mi = m 2 , as in Sec. II above [Eqs. (2.2a) and (2.2b)]. Even 
with unequal masses, at resonance, angle a is independent of the oscillator coupling k. 

With these results, Eqs. (3.6b) and (3.6c) give the normal-mode frequencies 



LU + = U) 

(4.4) 



«(mi + m 2 ) 



M 2 



w_ = Q = y u 
The ratio 7 = lj/Q can be written 

-('^)"' («) 

Here Q and 7 do depend on k. Since the masses must obey mi + m 2 > 2[i, with equality only for mi = m 2 , 
for given fi and u values frequency is greater, and ratio 7 is smaller, with unequal masses than in the 
symmetric case. 



B. Ground/coherent state: coefficient matrices. 

The behavior of the quantum state with oscillator initially in its unperturbed ground state, and 
oscillator #2 initially in a coherent state with classical amplitude x , follows immediately from the results 
of Sees. III.B-D above. With the frequencies appropriate to this case, Eq. (3.9e) gives a coefficient matrix S 
identical in form to that for the symmetric case, i.e., corresponding to wave function (2.13a-c). Matrices fii 



21 



and fi 2 are likewise the same as in the symmetric case; hence, so are U, U 1 , and V. The only matrices 
different in this case are 



f 



N= - : ( TOl m2 ] , (4.6a) 



its transpose iV T , its inverse 



1 



n- l = : ( M ™ 2 I , (4.6b) 



+m 2 ) TO i 
and the column matrix 

/ r / y/m 2 /m 1 sin(wi) \ 

2 = ^o, J" • (4.7) 

This follows from Eqs. (3.11c) and (4.3), expressed in terms of the original classical amplitude xq. 

The matrix products required are readily evaluated. They are 

, cos(u>t) — cos(fii) 
N- 1 U~ 1 VZ= m2X ° I , (4i 



mi + ra 2 



cos(ut) + (mi/mz) cos(Qt) 



N T UN 



/j(mi + m 2 )[cos 2 (fit) + 7 2 sin 2 (fit)] 

' ^ + m? [cos 2 (fit) + 7 2 sin 2 (fit)] -/i 2 (l - 7 2 ) sin 2 (fit) 

-M 2 (l - 7 2 ) sin 2 (fit) /i 2 + to 2 [cos 2 (fit) + 7 2 sin 2 (fit)] 



(4.9) 



n t \u + (v - n 2 )t7 _1 y]z = 



wsin(c<jt) — fisin(fit) 
"" J " ! ~ ' (m 2 /mi)u/sin(u;t) + fisin(fit) 



(4.10) 



and 



n t [u+(v - riiW^iv - n 2 )]jv = 



/z(mi + m 2 ) 



m( + y 2 [ cos 2 (fit) + — sin 2 (fit) 



1 /i 2 sin 2 (fit) 



1 



- 1 W sin 2 (fit) 



m\ + y 2 ^cos 2 (fit) + \ sin 2 (fit) 



/ 



(4.11) 
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C. Expectation values and uncertainties. 

The expectation values and uncertainties for position and momentum of oscillator #1 in this state are 
obtained from the above matrices via Eqs. (3.22a), (3.22b), (3.23), and (3.29). They yield 

(yi) = x a — — [cos(ujt) — cos(fii)] , (4.12a) 

mi + m 2 



h mi 



1 + _^[ cos 2(n t ) + 7 2 sin 2 (nt)} , (4.12b) 



y V 2m\uj m\ + m 2 \ mi 

(pi) = — mix — [wsin(wt) — fisin(fii)] , (4.12c) 

mi + m 2 



and = . [l + 2* f ^(flt) + J_ si n 2 (m)") 

p V 2 mi + m 2 L mi V 7 / 



(4.12d) 



The Gaussian probability distributions for position and momentum of oscillator #1 follow immediately from 
these results. 

The quantum behavior of oscillator #1 in this case resembles that in the in the symmetric case, but 
exhibits a broader range of possibilities. The momentum and position expectation values satisfy (pi) = 
mi d(yi)/dt. The position uncertainty oscillates between the unperturbed value and a smaller value: 



y 2miLd mi + m 2 y V 2miw 
The momentum uncertainty compensates: 



Jhrmu < (1) < % mi u mi +7 2 m 2 
V 2 ~ p ~ V 2 mi+m 2 ' 



The uncertainty product oscillates between the Heisenberg minimum value and a larger value- 



(4.15a) 



implying 



for mi > m 2 , or 



S<W<5a A + (4.15b) 



2 ~ v p ~ 2 y (mi +m 2 ) 2 V7 



^<4 1 M 1) <l(^ + 7) (4.15c) 

for mi < m 2 — always satisfying the Heisenberg Uncertainty Principle. However, unlike the symmetric case, 
for unequal-mass oscillators at resonance the "apparent loss of quantum nature," i.e., the quantum squeezing 
of the uncertainties, can be extreme. For a system with mi <C m 2 and strong coupling, i.e., 7 < 1, the 
minimum value of can be arbitrarily close to zero, the corresponding maximum value of arbitrarily 
large. Of course the squeezing is oscillatory, with and returning to their unperturbed quantum 
values y/h/ (2miw) and ymmiw/2 in each cycle. 
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V. CONCLUSIONS 



By virtue of its exact solubility, both classically and quantum-mechanically, a pair of linearly coupled 
harmonic oscillators proves to be a useful toy model for probing the interaction of a quantum system with its 
environment. Even in its ground state, this system displays effects of the coupling on the quantum nature of 
an oscillator: reduction of its position uncertainty below the unperturbed quantum value; a compensating 
increase in its momentum uncertainty, yielding an increase in the uncertainty product; breaking of the 
symmetry between position and momentum variables, a consequence of the position-dependent coupling. 
These cannot be characterized as finite-temperature effects. 

The interaction of a quantum oscillator with a classical one can be simulated via the coupled-oscillator 
pair in a quantum state with one oscillator initially in its unperturbed ground state, the other initially in a 
coherent or Glauber state incorporating classical behavior. The subsequent evolution of the wave function is 
calculated exactly, using ordinary harmonic-oscillator propagators for the normal modes of the system. The 
reduced probability distribution for the position of the initially quantum oscillator — a Gaussian distribution 
with time-dependent expectation value and uncertainity — and its position and momentum expectation values 
and uncertainties all follow from this wave function. The expectation values can be characterized as a "beat" 
amplitude between the normal modes. The behavior of the uncertainties, i.e., the quantum character of the 
oscillator, can be quite complicated. For oscillators with equal unperturbed frequencies, e.g., at resonance, 
this behavior can be described as a time-dependent quantum squeezing: The position uncertainty oscillates 
through values below the unperturbed value; depending on the relative masses of the original oscillators and 
the strength of the coupling, its lower bound can be arbitrarily close to zero. The momentum uncertainly 
oscillates through values larger than its unperturbed value. The uncertainty product oscillates through values 
larger than the Heisenberg minimum value. Naturally, the system never actually violates the Uncertainty 
Principle. 

Nonetheless, the behavior of the system is certainly reminiscent of the longstanding suggestion that 
a quantum system coupled to a classical one would "radiate away" its quantum nature. In detail, though, 
it more closely resembles the quantum squeezing encountered, e.g., in quantum fields coupled to classical 
backgrounds [8]. Although not explored in the present calculations, more involved effects connecting the 
quantum and classical worlds, such as decoherence, can be studied exactly and analytically in the two- 
oscillator system. 
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